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Abstract — This paper investigates the state feedback control problem of nonlinear networked control systems in 
presence of data packet dropout. The random data packet dropout is modeled as a Bernoulli distributed sequence 
with a known conditional probability distribution. The closed-loop system model is obtained through the state 
augmentation method. The sufficient conditions for the existence of a state feedback controller are derived by means 
of Lyapunov stability theory such that the closed-loop control system is exponentially stable in the mean-square sense. 

The maximum bound of the nonlinear item of the control system is also obtained. Then a linear matrix inequality 
(LMI) based approach for designing such a controller which can be solved by Matlab LMI tool box is presented. 
Finally, a numerical simulation example is used to demonstrate the effectiveness of the proposed method. 
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I. Introduction 

In networked control system (NCS), control loops are closed by real-time networks. Such kind of control systems bring a lot 
of advantages, such as reduced system wiring, simple system maintenance and low cost. In the field of industrial control, 
remote medical treatment, military and other fields, NCS has been widely used [1-2]. However, the application of 
communication networks in feedback control loops makes the NCS analysis and synthesis complex because data packets 
through networks suffer not only transmission delays but also data packet drop out [3-4] .Many scholars have investigated the 
effect of data packet dropout on the performance of the control system. 

In recent years, there are many results for the controller design of linear network control systems respect to the data packet 
dropout. For the NCS with random data packet drop out, literature [5-6] modeled NCS as a kind of Markov jump system, and 
the mode dependent controller is designed. In literature [7], if the data packet drop out accrues, the input of the actuator was 
considered to be zero, and the optimal estimation of NCS was studied. For a class of singular controlled plant, literature [8] 
modeled closed-loop system as a kind of asynchronous dynamical system, and the H m control problem was investigated. 
However, there is not much work reported about nonlinear NCS. The state feedback control problem for nonlinear NCS is 
still very limited. 

In this paper, the state feedback control problem for a kind of nonlinear NCS is investigated. Assuming that the data packet 
dropout obeys Bernoulli distribution, the closed-loop system is modeled as a class of control system which contains 
stochastic variable. The sufficient condition and design method of the state feedback controller for a class of nonlinear NCS 
with random data packet drop out are given. Simulation results show that the method proposed in this paper is effective. 

II. Problem Formulation 

A nonlinear controlled plant is given by: 

x(k + 1) = Ax(k) + Bulk) + / ( k , x(k)) ( 1 ) 

Where x(k ) e R" is system state vector and ll(k ) e R"" is the system control input. A and B are all real constant matrices. 

/ (k,x(k)) is nonlinear disturbance, satisfying the following Lipschitz condition: 

f T (k,x(k))f(k,x(k)) < S 2 x T (k)H r Hx{k) ( 2 ) 

Where His real constant matrix, ando > Ois the bound of the nonlinear disturbance. 

We consider the circumstance that data packet drop out occurs between sensor and controller. A stochastic variable is defined 
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as follows to describe the process of data packet drop out : 

[ 0 data packet dropout occurs 
[1 otherwise 

The mathematical expectation and variance of a(k) are as follows: 

Prob{«(k) = 1} = E{a(k)} := a , Prob{a(k) = 0} = 1 — E[a(k)} := 1 —a , Var{o'(k)} = E{(a(k)-a) 2 } = a{l-a) := b 2 . 

When data packet dropout occurs, the actuator continues to use the previous cycle control signal, so NCS with data packet 
dropout can be described as: 

x(k + 1) = Ax(k) + a(k)Bu(k) + (1 - a(k))Bu(k - 1) + f(k, x(k)) (3) 

Assume the state feedback control law is: 

u(k) = Kx(k) (4) 


Substituting (4) into (3), we have the following closed -loop system: 

x(k + 1) = (A + a(k)BK)x(k) + (1 - a(k))BKx(k - 1) + f(k,x{k)) 

~lT 


(5) 


written as: 


At sampling timek , if we augment the state-variable as %(k) T =[ ~x(k) T x(k- l) r ] .the closed-loop system (5) can be 

(6) 

, e = i j , F(k, m) = f(.k, cm) ,g=[i o] . 


Where 


A (k) = 


%(k + 1) = A (k)£(k) + EF(k, %{k)) 
1 ~ 


A + a(k)B (1 -a{k))B 
I 0 


Before proceeding, we need the following definition. 

Definition 1 The closed-loop NCS (6) is said to be exponentially mean-square stable if there exist constants /U > 0 and 
0 < p < 1 such that s{ II II 2 } < fup k s{ II E, k II 2 } holds for all. 

Lemma 1 [9] LetTJ, T 2 be symmetric matrices. The conditions on X T T x X < Osuch that x' T n X < 0 hold if there exist a scalar 
£ > 0 such that 7J — sT 2 < 0 . 


III. Controller Design 


Theorem 1 Take as given communication channel parameters Prob{«(k) = 0} = l — a .The closed-loop NCS(6) is 
exponentially mean-square stable if there exist positive definite matrices P,Q and nonnegative real scalar £ > Osuch that the 
following matrix inequality holds: 


E n +eS 2 H T H 

“21 

“31 


22 

"'32 


* 

* 

, ~el 


<0 


(7) 


where 

S„ = (BK) t P(A + aBK) + b 2 ( BK) T PBK + Q - P , S 21 = (1 - a)( BK) r P(A + aBK) - b 2 ( BK) r PBK , S 31 = IP(A + aBK) , 


E 22 =(l-a) 2 (BK) T PBK+b 2 (BK) T PBK-Q ,Z 32 =(l-a)IPBK ,S 32 = P . 
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Proof: 

Define a Lyapunov function 

v<jc)=z(k) T sm=t(ky 


P 0 
0 Q 


m T = 


x(k) T Px(k ) + x(k - 1) T Qx(k - 1) 


Where P and Q are positive definite matrices, we have: 


E{V k+1 }-E{V k } 

— E{ [( A + a(k)BK)x(k) + (1 - a(k))BKx(k - 1) + f(k, x(k)] T P 
[( A + a(k)BK)x(k) + (1 - a(k))BKx(k - 1) + f(k, x(/c)] + 
x(k) T Qx(k) - x(k) T Px(k ) - x(k - l) r Qx(k - 1) 

= E{ [(A + aBK)x(k) + (a(k) -a)BKx(k) -(a(k) -a)BKx(k -1) + (l-a)BKx(k — 1) + f(k,x(k)] T P 
[( A + aBK)x(k) + ( a(k ) — a)BKx(k) - ( a(k ) - a)BKx(k — 1) + (1 -a)BKx(k - 1) + f(k, x(/t)] 
x(k) T Qx(k) — x(k) T Px(k) — x(k -l) r Qx(k -1) 


=m r 

Where </>{k) T = [x(^) r x(k-l) T f(k,x(k)) T ~^ 
From (2), we can get: 

m 1 


^li 

"'21 “22 
■"31 “32 


mk) 


-8 PI PI * * 
0 0 * 
0 0 1 


</>(k) < 0 


( 8 ) 


(9) 


According to lemma 1, if there exists a scalar e > 0 such that 


E n +s8 2 H T H 

“21 

“31 


■"22 

"'32 


* 

* 

, -el 


<0 


then the following inequality subject to (9) holds. 



( 10 ) 


Since S > 0 ,it is clear that 


Kin(s)£{\\ri k \\ 2 } < s{ri T k Sr] k } < A max (sMH^II 2 }. Because e{f k+ ^^ k+l }-£{^ T k S^ k } < 0 , there existsO < y < A max (S)such 
that e{f[ +1 S^ k+1 ] < -YlUk II 2 } + ^k S ^k) ^ (l - t Y m ) Consequently, we can get 

' A max WJ' 


A mm (sM\\ri k \\ 2 }<s{riSB k }<(l- 


r 

AS) 


M £sz k }< 
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(i - -7-^)* } < 4™ (S)( i - t-^)* 11 4 11 2 1 


KJS) 


Kr^S) 


That is £{\\ 4 \\ 2 }</jp k s{\\^ k II 2 }, 


where p = ^ max ^ , p = 1 — — , 

^(S) ^(S) 


Theorem 2 Take as given communication channel parameters Prob{«(k) = 0} = 1 - a , the systems (6) is exponentially mean 
square stable if there exist positive definite matrices M>0,N>0, matrix Y and a scalar (p > 0 satisfying the following LMI 

'£ *‘ 


n a 


<0 


where 


Z = diag(N -M ,-N,-I ) , A = diag(-M ,—M ,-tpI ) , fl = 

Proof: 

(7) Is equivalent to 


AM (1 - a) BY I 
bBY -bBY 0 

HM 0 0 


( 11 ) 


and the controller is given by u(k) = YM ' x(k ) . 


i n +8 H H * * 

£7 * 


J 21 

■ J 31 


, -I 


<0 


( 12 ) 


Where 


S n = (BK) r P(A + aBK ) + b 2 ( BK) T PBK +Q-P, S 21 = (1 - a)(BK) T P(A + ciBK) - b 2 (BK) T PBK , 
S 31 = 1P(A + aBK) , S 22 = (1 - a) 2 (BK) T PBK + b 2 ( BK) T PBK -Q, S 32 = (1 - a)IPBK , S 32 = P , 

P= P -,Q = ± 

£ £ 


And (12) can be written as 


b 2 ( BK) T PBK + Q-P + 8 2 H T H 

* * 


~ (A + aBK) 1 " 

-b 2 (BK) T PBK 

b 2 (BK) T PBK-Q * 

+ 

(1 -a)(BK) T 

IP( A + aBK) 

0 -I 


I 


P[A + aBK (1 -a)(BK) /]<0. 


By Schur complement, one can obtain 


b 2 (BK) T PBK + Q-P + S 2 H T H 
-b 2 (BK) T PBK 
IP( A + aBK) 

A + aBK 


b (BK) PBK- Q * * 

0 -/ * 
(1 -a)BK I -F 


<0 


Similarly, using Schur again we can get 


A * 
^11 

A 2 i A 22 


<0 


(13) 


Page | 39 
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Where 


A n = diag(Q- P.-Q.-I) , A 22 = dicig(-P,-P,-<pI) , A 21 


* * 

bBK -bBK * . 
H 0 0 


Performing a congruence transformation to (13) by diag{P 1 , P *,/,/, /,/,/}, and let M =P 1 , N = MQM ,Y = KM 


5 2 = (p , (1 1) can be obtained, which complete the proof. 


Corollary 1 If the following optimization problem has solution 

min (p s. t. (11) 

M>o,w>o,<p>o,r 

and the solution is (p min ,the maximum bound of nonlinear disturbance in (2)is 8 max = yj(p min . 

Remark 1 the above convex optimization with LMI constraints can be solved by Matlab LMI toolbox objective function 
minimization problem solver mincx. 


IV. Simulation Example 


Consider the following nonlinear controlled plant 


x(k + 1) = 


1 

0.3 


-2 

0.6 


x(k ) - 


-1 

0.2 


u(k) - 


0. lx, (k)sin x t (k) 
0.2x 2 ( k ) cos x, ( k ) 


Assume that the initial state of the system is x(0) = [3 -l] r and the Probability of the data packet dropout is 
Prob{«(k) = 0} = 0.2 and the matrix H is a unit matrix. According to Theorem 2, and using Matlab LMI tool, we can get 


Y = [0.4610 -0.2155], M 


40.0981 

29.5003 


29.5003 

63.9610 


. Hence, the controller gain matrix is 


K = FM 1 = [0.0212 -0.0131] . 


By Matlab LMI toolbox objective function minimization problem solver mincx, the maximum bound of nonlinear 
disturbance in (2) is obtained as 8^ = 0.2277 . The process of data packet drop out and the state trajectory of the closed-loop 
system is shown by the Fig. 1 and Fig. 2, respectively. 



Fig.1 The random data packet dropout distribution 
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Fig.2 The state trajectory of the closed-loop system 

From the figure above, we can see that the closed-loop system is still stable under the designed controller even if there exists 
data packet dropout. 

V. Conclusion 

In this paper, the sufficient condition and design method of the state feedback controller for a class of nonlinear NCS with 
random data packet drop out are given. The controller gain matrix and the maximum bound of the nonlinear item of the 
control system can be obtained by use of the relevant tool box of Matlab. Simulation results show that the method proposed 
in this paper is effective. 
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